An Improved Experimental Limit on the Electric-Dipole Moment of the Neutron 
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An experimental search for an electric-dipole moment (EDM) of the neutron has been carried out 
at the Institut Laue-Langevin (ILL), Grenoble. Spurious signals from magnetic-field fluctuations 
were reduced to insignificance by the use of a cohabiting atomic-mercury magnetometer. Systematic 
uncertainties, including geometric-phase-induced false EDMs, have been carefully studied. Two 
independent approaches to the analysis have been adopted. The overall results may be interpreted 
as an upper limit on the absolute value of the neutron EDM of \dn\ < 2.9 x 10~^^ ecm (90% CL). 

PACS numbers; 13.40.Em, 07.55. Ge, ll.30.Er, 14.20.Dh 



I. INTRODUCTION 

Measurements of particle electric-dipole moments 
(EDMs) P, 0, 13 are of significant interest because they 
provide some of the tightest constraints on extensions to 
the Standard Model, such as supersymmetry, that at- 
tempt to explain the mechanisms underlying CP viola- 

tion iiiaiiiiaini 

This neutron-EDM experiment, and the performance 
of its cohabiting mercury magnetometer, have been dis- 
cussed in earlier publications The final result pre- 
sented in this Letter incorporates a comprehensive analy- 
sis of systematic errors, some of which were undiscovered 
at the time of the earlier measurements. 



II. EDM MEASUREMENT TECHNIQUE 

The measurement was made with ultracold neutrons 
(UCNs) stored in a trap permeated by uniform E- and B- 
fields. This adds terms — /x„ -B and — d„ -E to the Hamil- 
tonian determining the states of the neutron. Given par- 
allel E and B fields, the Larmor frequency with which 
the neutron spin polarization precesses about the field di- 
rection is 

= |2^„B + 2d„£;|. (1) 

For antiparallel fields, hi'^^ = |2/i,ii3 — 2d„i?|. Thus the 
experiment aimed to measure any shift in the transition 
frequency v as an applied E field alternated between be- 
ing parallel and then antiparallel to B. 

A schematic of the apparatus is shown in Figure ^ 
The UCNs were prepared in a spin-polarized state by 
transmission through a thin, magnetized iron foil, and 
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entered a cylindrical 21-litcr trap within a 1 /iT uniform 
vertical magnetic field Bq. 

Approximately 20 s were needed to fill the trap with 
neutrons, after which the entrance door was closed pneu- 
matically. The electric field, of approximately 10 kV/cm, 
was generated by applying high voltage to the electrode 
that constituted the roof of the trap, while keeping the 
fioor electrode grounded. The electrodes were made of 
diamond-like-carbon coated Al, and the side wall was 
SiOs. 

The transition frequency v of the neutrons was mea- 
sured using the Ramsey separated oscillatory field mag- 
netic resonance method. During the storage period, the 
neutrons interacted coherently with two 2 s intervals of 
oscillating magnetic field having a chosen frequency close 
to the Larmor frequency. The two intervals were sepa- 
rated by a period T = 130 s of free precession. The last 
step was to count the number of neutrons 7V| and iV^ 
that finished in each of the two polarization states. This 
was achieved by opening the entrance door to the trap 
and allowing the neutrons to fall down onto the polariz- 
ing foil, which then acted as a spin analyzer. Only those 
in the initial spin state could pass through to the detec- 
tor, which was a proportional counter in which neutrons 
were detected via the reaction n+'^He ^■^H-f-p. During 
one half of the counting period, an r.f. magnetic field 
was applied in the region above the polarizing foil; this 
fiipped the spins of the neutrons, thereby also allowing 
those in the opposite spin state to be counted. Each 
batch cycle yielded about 14,000 UCN counts. Within a 
run the data-taking operations were cycled continuously 
for 1-2 days. Periodically, after a preset number (nor- 
mally 16) of batches, the direction of E was reversed. 
All other settings were held constant during a run. Ev- 
ery 10-20 runs, Bq was reversed so that half of the full 
data set was taken with Bq upwards and half with Bq 
downwards. We adopt a system as in ^T§\ where the k 
vector of our z axis follows the direction of Bq. Hence, Bq 
is always positive, while the gravitational displacement of 
the UCNs changes sign. 
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FIG. 1: (Color online) Experimental apparatus 



The magnetometer used the precession frequency of 
/ = 1/2 atoms of ^^^Hg (3 x 10^*' atoms/cm^; firi/f^Hg = 
In/lHg = —3.842) stored simultaneously in the same 
trap as the neutrons. Using Eq. for both UCNs and 
Hg, and assuming that both experience the same B, we 
find that to first order in the EDMs d, 
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For each run, d„ 



(2) 

was obtained from a linear fit to the 
ratio Vn/vHg versus E. Eq. shows that dmeas contains 

has been shown 
, so it introduces 



a contribution from dng- The true d 
to be (-1.06±0.49±0.40) x IQ-^s e cm 
a systematic error of (— 0.4±0.3) x 10"^"^ ecm into d; 



(warmer) Hg atoms, so an observed shift of Vn/^Hg away 
from \^n/lHg\ gives a measure of the volume-averaged 
{dBoJdz)v, via the resuh (Eq. 86 in fl^) 
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where the -I- sign corresponds to Bq downwards. 

In this experiment the contribution of the GP effect in 
the Hg to dmeas IS 50 times larger than the GP effect of 
the UCNs. Writing the GP false contribution to dmeas 
from the Hg as dn_Hg.f, it is shown in (Eq. 87) that 



±- 



HnlHgl 



Ahc'^ 



iRa~l) = ±fc-(i?,-l), (4) 



where rg is the trap radius and the -I- sign again corre- 
sponds to Bq downwards. It follows that we can write 



= d' 



''-nMgJ 



^dl,±k-{B.a-Rao), (5) 



where is the true d^ plus all other systematic ef- 
fects discussed below, and Rao is the value of Ra where 
dBz/dz = 0. Eq.O defines two straight lines, one with 
positive slope for Bq down and one with a negative slope 
for Bq up. The crossing point (i?aO,dJi) provides an es- 
timator of d'^ free of dn,Hg,f- 

Each run was made at a chosen value of Ra by pre- 
adjusting currents in field-trimming coils. Figure[21shows 
the data (binned for clarity) for dmeas as a function of 
Ra for each direction of Bq. The lines represent a least- 
squares fit to all 554 of the (unbinned) run results, using 
as free parameters the two intercepts and a common ab- 
solute slope k. This yields x^/^ — 652/551 and k — 
(1.90 ± 0.25) x 10-26 ecm/ppm, which is within 1.3cr 
of the expected value of (1.57 ± 0.08) x 10-^^ ecm/ppm 
from Eq. Q . The slope k can be altered by a few percent 
(although still remaining highly symmetric under Bq re- 
versal) by various mechanisms including the UCNs' own 
GP signal (a 2% effect); uncertainty in Ah (4%); a slight 
reduction in mean free path due to cavities and grooves 
in the electrodes as well as to the presence of lO""^ torr of 
He gas to prevent sparks [T3L [l3| (1%); and asymmetric 
surface relaxation of the Hg (up to 5%). 



Anomalous GP effects 



III. GEOMETRIC-PHASE EFFECTS 

To the true d„ and dng within dmeas there will also 
be added coefficients of fractional shifts in and VHg, 
from other causes, which are linear in E and thus consti- 
tute additional systematic errors. The most important 
of these involves a geometric phase (GP) arising when 
the trapped particles experience a gradient dBg^/dz in 
the presence of E [l^ ■ Fortunately, the centre of gravity 
of our UCNs is Ah = 0.28 cm lower than that of the 



There are some processes that can interfere with the 
above GP error removal ~ essentially any process that 
changes Ra and/or dn^Hg.f without conforming to the 
ratio between the two given by Eq. (0J|, and where, in 
addition, the changes differ with the direction of Bq. 

First, there are several processes that shift Ra but 
not dn,Hg,f- Changing ^n/lHg, for example, shifts the 
two lines of Figure 2 in the same direction by the 
same amount, leaving d'^ unaffected but changing iJ^o- 
We note at this point that our final Rao is consistent 
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FIG. 2: (Color online) Measured EDM (binned data) as a 
function of the relative frequency shift of neutrons and Hg. 



Second, there are processes that generate an enhanced 
dn,HgJ- The field of a permanent magnetic dipole 
(PMD) close to the trap makes a non-uniform dBz/dz 
and adds an (enhanced) GP ddip To'l to dn^Hgj, but shifts 
Ra in accord with Eq. ((SJ and in opposite senses for the 
two Bq directions. The two changes are in a ratio greater 
than that given by Eq. This shifts the lines of Figure 
2 upwards, adding {ddip — d^) to dj^, where d^ is the pre- 
diction of Eq. Q for the {dBz/dz)v of this PMD. Our 
fiuxgate magnetometer surveys of the trap cannot rule 
out PMD fields of less than 1 nT at 2 cm from the inner 
surface. Large areas of the trap are Si02 or Al, backed 
by large voids, and do not come under suspicion; but the 
Hg and UCN doors involve a heterogeneous collection 
of small parts close to the trap. We allow a d'^ uncer- 
tainty of ±6.0 X 10~^^ e cm to allow for an undetected 1 
nT PMD at the Hg door. In the case of the UCN door 
we have better diagnostics. The trap used when taking 
EDM data has a small cavity in the lower electrode, 4.0 
cm deep and 6.8 cm in diameter, sealed from below by 
the door. A PMD in the door mechanism can contribute 
a dBz/dz field to the cavity and to the rest of the trap. 
The field in the cavity contributes a strong shift in Ra, 
while contributing negligibly via the GP to dn,Hg,f due 
to the small cavity radius. As a result, the lines of Fig- 
ure 2 are shifted in opposite directions, again adding a 
systematic error to d'^. 



B. Auxiliary measurements 



with the "fn/ynn value from the literature (Icr = 1.7 
ppm) [T5I llq| . after allowing for our observed Bj; and 
By fields with finite dB^/dy and/or dBy/dx but with 
(dBx/dx + dBy/dy) = = —dBz/dz (e.g., a quadrupole 
aligned with z, with B^, — qy,By = qx), which cause 
Ra to increase without contributing to dn,Hgj- Any 
change in such fields when Bq is reversed can result in 
a differential shift of all Ra values, and thus of the two 
lines, thereby changing Rao and d'„. Below, we describe 
our measurements of the differential shift in Rao and state 
the correction to d'^. 

B-field averaging in the trap is affected by localised 
loss of UCN and Hg particles, and by polarization loss in 
the presence of the 10""^ fractional Bq inhomogeneities, 
which may change with Bo direction. However, we es- 
timate that the resulting Ra shifts are < 0.1 ppm and 
< 0.01 ppm for the UCN and Hg respectively, and that 
they will be indistinguishable from the quadrupole shifts. 

Light shifts 17, 18] in h'Hg will shift Ra- They are 
produced by any small component, parallel to Bq, of the 
204jjg ppQ^jg light beam passing through the precessing 
^^^Hg atoms. This component, and the Ra shift, reverse 
sign on reversal of Bq. A slight dependence of Ra on the 
incident light intensity was indeed found, the magnitude 
^ 0.2 ppm being in agreement with theory. A correction 
to dmeas was made on a run-by-run basis, leading to an 
overall correction of (3.5 ± 0.8) x 10~^^ ecm. 



Our additional diagnostics came from separate Vn, t^Hg 
and Ra measurements (without E fields) in two auxil- 
iary traps having roofs that could be raised or lowered to 
change the height H. The traps were built on the same 
lower electrode and door mechanism that were used for 
EDM data taking. Assuming Bz{z) — bo + biz + b2z'^ , one 
can show that {dBz/dz)v = for a trap height H when 
VHg is the same for roof settings at H/2 and H . With 
such a field established, Rao was measured for several 
heights H. The resulting forms of Rao{H) for Bq up and 
Bq down led to the conclusion that there was a dipole 
field of strength I nT penetrating into the door cavity. 
They also confirmed the presence of quadrupole fields, 
which the fluxgate scans showed vary little with H. Po- 
larization data from the EDM runs further substantiated 
this: UCNs of different energies have different Ahs, so 
the surviving UCN polarization decreases in proportion 
to (dBzdz)^, and the values of Ran where the Bq up and 
down polarizations peak are in excellent agreement with 
the auxiliary trap results. 

The first auxiliary trap used had a smaller radius, 18.5 
cm rather than 23.5 cm, and a cavity depth of 6.0 rather 
than 4.0 cm. These differences amplify the Ra shifts 
from a dipole field in the cavity by 1.50 and reduce the 
quadrupole shifts by a factor 1/1.8. Our systematic error 
correction to to allow for the combined door dipole 
and quadrupole fields is (0.69 ± 0.28) x lO'^e gem. The 
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Effect 


Shift 


a 


Door cavity dipole 


-5.6 


2.0 


Other dipole fields 


0.0 


6.0 


Quadrupole difference 


-1.3 


2.0 


V X E translational 


0.0 


0.03 


V X E rotational 


0.0 


1.0 


Second-order v x E 


0.0 


0.02 


UHg light shift (geo phase) 


3.5 


0.8 


VHg light shift (direct) 


0.0 


0.2 


Uncompensated B drift 


0.0 


2.4 


Hg atom EDM 


-0.4 


0.3 


Electric forces 


0.0 


0.4 


Leakage currents 


0.0 


0.1 


AC fields 


0.0 


0.01 


Total 


-3.8 


7.2 



TABLE I: Summary of systematic errors and their uncertain- 
ties, in units of lO^'^'^ ecm. 

auxiliary trap used to measure Ah by obtaining Ra as 
a (linear) function of a scries of known {dB^J dz)v was 
made as similar to the data-taking trap as possible in 
dimensions and materials so as to reproduce the same 
UCN velocity spectrum. 



IV. OTHER SYSTEMATIC ERRORS 

We now consider the systematic errors not involving 
GP. 



A. V X E effects 

If the UCN have a net translational motion, any per- 
pendicular component of the E field will be seen in their 
rest frame as a combination of E and B fields. During 
the 130 s Ramsey measurement period, the UCN ensem- 
ble may warm slightly due to vibrations. This causes the 
center of mass to rise by up to 0.1 cm. If the volume- 
averaged angles between E, B and v arc each as high as 
0.05 radians, the induced false EDM will be 0.03 x lO'^^ 
ecm. 

In a similar manner to the translational effect, any net 
rotational flow of the UCN in conjunction with a radial 
component of the E field may lead to an induced EDM 
signal. However, any such flow of UCN is expected to 
be attenuated by wall collisions before the first Ramsey 
pulse is applied. We calculate that the maximum error 
to be expected from this source and from higher-order 
V X E effects is below 1 x 10^-^^ e cm. 



B. Direct light shift 

Analysis of the data suggested a possible small corre- 
lation between the intensity of the Hg reading light and 



the value of the applied E field. Through the light shift, 
this could directly create an EDM signal. However, to 
within its uncertainty, the dependence of Ra on the light 
intensity has been removed. The residual systematic un- 
certainty from this source is 0.2 x 10"^'' ecm. 



C. Uncompensated magnetic field fluctuations 

There may also be residual effects from B field fluc- 
tuations. For example, a dipole-like field B^ originating 
from the /z-metal in the region of the HV feedthrough 
would be sensed by both neutrons and Hg, but with a 
difference given by SB^/Bd = 3Ah/r where r ~ 55 cm 
is the distance from the source of the field to the center 
of the bottle. Thus, fluctuations in B that are corre- 
lated with the HV can be expected to be compensated 
up to a factor of about 70. In order to study this, the 
Hg and neutron channels were analysed independently. 
The neutrons yielded an EDM signal of (17 ± 4) x lO^^e 
ecm; the Hg, once the GP contribution (as calculated 
from the average i?a — 1 at which the data were taken) 
was subtracted, yielded (-3.9 ±0.8) x 10~^^ ecm. These 
results are consistent with a common source of magnetic 
fluctuations correlated with the HV. We therefore expect 
the Hg compensation to shield us from this systematic 
effect to a level of 17 x 10-26/70 = 2.4 x lO'^^ ecm. 



D. Electric forces 

Another possible source of systematic error arises 
from electrostatic forces, which may move the electrodes 
slightly. In conjunction with a magnetic field gradient, an 
HV-dependent shift in the ratio would then appear. This 
was sought by looking for an EDM-like signal but with 
a frequency shift proportional to |E| instead of E. The 
|E| signal was consistent with zero, with an uncertainty 
of 4 x 10~^^ e cm. If the HV magnitudes were slightly 
different for the two signs of E, this effect would generate 
a false EDM signal. Study of the measured HV and of 
the charging currents show that the HV magnitude was 
the same for both polarities to within about 1%. This 
systematic uncertainty is therefore 0.4 x 10"^^ e cm. 



E. Leakage currents and sparks 

Analysis of the EDM as a function of leakage current 
shows no measurable effect. Leakage currents are typ- 
ically of order 1 nA. If this current were to travel 10 
cm azimuthally aroimd the bottle, the resultant B field 
would result in an apparent EDM of 0.1 x 10"^^ ecm. 
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F. HV AC ripple 

AC fields are another possible cause of concern. There 
is no differential ripple visible on the HV at the level of 
a few volts. Sampling is done at 5 Hz with a bandwidth 
of 20 kHz, so any 50 Hz ripple would show up as beats. 
This is certainly absent at the level of 50 V, which would 
give a false EDM of 0.01 x 10"" g cm. 

Low-frequency AC fields were sought by means of a 
pickup coil in conjunction with a phase-sensitive detec- 
tor. Shifts in Ra from this source at the level of 0.02 
ppm could not be ruled out. Cancellations in the corre- 
sponding EDM signal from reversals of the electric and 
magnetic fields would reduce any net contribution to be- 
low the level of 0.01 x 10^^'' ecm. 



V. ANALYSIS AND RESULTS 

Two approaches were adopted in studying the data. 
In the first, more straightforward, analysis, only the 293 
runs with an uncertainty on Ra of less than 0.05 and 
with -10 ppm < {Ra — 1) < 21 ppm were accepted. As 
the two resulting average -Ra — 1 values (8.948, 8.943 ppm 
for Bq up, down respectively) were almost identical, with 
approximately equal amounts of data in each field direc- 
tion, a simple weighted average of the dmeas data was 
used as an estimator of (i„. The value obtained was 

dn = (-0.2 ± 1.6 (stat)) x 10"^^ e cm {x^/iy = 1.24), 

with an additional uniformly distributed systematic un- 
certainty of ±1.2 X 10~^^ ecm allocated to it in order 
to accommodate any potential systematic biases arising 
from the effects listed in Table 1. The resulting dis- 
tribution of possible values implies an upper limit of 



|d„| < 2.9 X 10-26 ecm (90% CL). 

The second analysis began with the Figure El fitted- 
lines crossing-point value = (—0.55 ± 1.51) x lO-^^ 
e cm (which includes the run-by-run light shift correc- 
tion) and then applied the systematic-error corrections 
given in Table 1. The final result from this approach is 

dn (+0.2 ± 1.5 (stat) ± 0.7 (syst)) x 10^^^ ecm. 

In this analysis, the systematic uncertainty is normally 
distributed. The result implies an upper limit of \dn\ < 
2.8 X 10^26 g(,j^ (90% CL)^ 



VI. CONCLUSION 

The data set analysed here, which excludes data that 
have already been published [l| , incorporates all neutron 
EDM measurements undertaken between the autumn of 
1998 and the end of 2002. The results overall may be 
interpreted as an upper limit on the absolute value of 
the neutron EDM of |d„| < 2.9 x IQ-^e ecm (90% CL). 
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